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Abstract 

We reinterpret the generic CDF charged massive particle limit to obtain a limit on 
the mass of a stable or long-lived gluino. Various sources of uncertainty are examined. 
The i?-hadron spectrum and scattering cross sections are modeled based on known 
low-energy hadron physics and the resultant uncertainties are quantified and found 
to be small compared to uncertainties from the scale dependence of the NLO pQCD 
production cross sections. The largest uncertainty in the limit comes from the unknown 
squark mass: when the squark - gluino mass splitting is small, we obtain a gluino mass 
limit of 407 GeV, while in the limit of heavy squarks the gluino mass limit is 397 GeV. 
For arbitrary (degenerate) squark masses, we obtain a lower limit of 322 GeV on the 
gluino mass. These limits apply for any gluino lifetime longer than ~ 30 ns, and are the 
most stringent limits for such a long-lived or stable gluino. 



1 Introduction 

The observation of exotic stable massive particles (SMPs) which can be detected by their 
interactions in a detector would be of fundamental significance. SMPs are features of a 
number of scenarios of physics beyond the Standard Model, such as theories of supersymmetry 
(SUSY) and extra dimensions [TJ. Searches have therefore been carried out at colliders, in 
cosmic rays and matter [TJ [5] . Collider searches for coloured SMPs - i?-hadrons in the context 
of SUSY - present an additional challenge in interpreting the experimental observations, 
owing to uncertainties in the mass spectra and scattering of the colour-singlet i?-hadrons. 
This is of particular concern for i?-hadrons consisting of a heavy colour octet combined with a 
gluon or qq in a colour-octet state, since these systems have no direct analog among Standard 
Model hadrons [3]. Consequently, there have been no dedicated searches for stable gluino 
i?-hadrons since the LEP era [U[5]. In this paper, we outline an approach to modeling the 
mass spectrum and the scattering of gluino-based i?-hadrons. This provides the theoretical 
and phenomenological foundation for a new method to search for long-lived or stable massive 
gluinos. Using our method, we reinterpret a recent CDF limit [5] to give a new lower limit 
on the mass of a stable or long-lived gluino. 

There have been dedicated searches for unstable massive gluinos since the LEP era. D0 per- 
formed a search for unstable gluinos that stop in the detector and decay [7J. While the 
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present paper was in review, CMS presented results of an analogous search at the LHC [8] . 
Searches for decaying gluinos are highly dependent on the lifetime of the gluino, whereas our 
approach only requires that the gluino have a lifetime long enough that it traverses the de- 
tector without decaying in flight, r > 30 ns, a condition satisfied by all the models considered 
in stopped gluino searches. Thus our limits can be directly compared to those studies but 
have a larger range of applicability. For brevity below, we call a particle stable if it does not 
decay during its traversal of a detector. 

There are two main approaches to search for non-decaying strongly interacting SMPs at a 
collidefl One strategy exploits the expected anomalous energy loss of a SMP as it propagates 
through an inner tracking system next to the beam-line, while the other relies on the SMP- 
speed (,8 <C 1) and possible penetrating behaviour whereby a SMP can be measured as a 
slow-moving object in an outer muon tracking system. Typically, experiments at LEP [H [5] 
and HERA [9 exploited the former, while the Tevatron experiments |6j HQ] relied on the latter 
signature. However, since i?-hadrons interact hadronically in dense calorimeter material 
positioned between the aforementioned tracking systems [TTJ [TJ1 [T3] , the i?-hadron charge 
can change between the inner detector and when it reaches the outer muon system. The 
charge of the i?-hadron in the outer muon system depends both on the scattering mechanism 
and the i?-hadron mass hierarchies. A comprehensive search for i?-hadrons at a collider 
must consider i?-hadron models in which the i?-hadron i) is dominantly neutral throughout 
its passage in a detector, ii) can be charged in the inner system but is dominantly neutral in 
the outer tracking systems, and Hi) can be charged in both the inner and outer trackers. Since 
the current strongly-interacting SMP-limits from the Tevatron assume the latter scenario, we 
provide here a gluino mass limit applicable to models which predict the muon-like signature 
Hi). 

This paper is organised as follows. First we outline expectations for the gluino- containing 
i?-hadron mass spectra. This is followed by a description of the simulation of i?-hadron 
scattering in matter. Using the material composition of the CDF detector and fixed order 
QCD-models, the proportions of gluino i?-hadrons which would have passed the CDF SMP- 
selections are then estimated and mass limits extracted. 

2 Spectrum of R-hadron masses 

To consider the experimental signatures of i?-hadrons in a detector we need to postulate a 
mass spectrum for the i?-mesons and the i?-baryons. The masses of the i?-hadrons have been 
calculated using the MIT bag model [18] for i?-mesons pjJ] and i?-baryons [2oJl. The mass 
splittings within a given i?-hadron multiplet are governed by the same QCD interactions that 
are responsible for the mass splittings of ordinary mesons and baryons. These are surprisingly 
well-described by an effective one-gluon-exchange interaction, and many unknown parameters 
for the i?-hadron mass splittings can be fixed by splittings in the analogous ordinary hadrons. 
To the extent that parameters are determined phenomenologically, splittings should be fairly 
model- independent. 

An important feature of the i?-baryon mass spectrum, first pointed out by Farrar [3], is 
that the lightest i?-baryon is the neutral flavor-singlet udsg, denoted S°. Although the 
constituent quarks are the same, this is a distinct state from the flavor-octet A . Due to the 

1 We only consider strongly interacting SMPs in the remainder of this work, and denote them hereafter as 
simply SMPs for brevity. 

2 Latticc calculations have been performed for R-mesons and gluino-balls, e.g., 1211 . but the bag-model 
results are consistent with the lattice calculations and easier to use. 
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particularly strong hypcrfine attraction in the flavor-singlet channel, the S° is lighter than the 
i?-proton. There is no flavor-singlet ordinary baryon analog of the S , due to Fermi statistics 
when the quarks form a colour-singlet. As noted above, the QCD hyperfine interactions that 
govern the R-baryon spectrum are well-constrained by the spectrum of ordinary hadrons 
and our knowledge of low-energy QCD, making it highly unlikely that the S° would not be 
the lightest state. The signal in our analysis would be bigger for a given gluino mass if the 
lightest R-baryon were charged, so doing the analysis with the most plausible mass spectrum 
(the S° being lightest) gives the most conservative gluino mass limits. 

The i?-meson spectrum is important for our analysis, because the initial i?-hadrons resulting 
from hadronization of a gluino are predominantly i?-mesons. i?-mesons are formed by the 
combination of the gluino with a colour octet qq pair. We must also consider the gluino-ball 
gg on a similar footing as it mixes directly with the flavor singlet _R-meson states. The states 
can be broken down by first considering the qq state and then combining this with the spin-^ 
g. The combination of the J c = 0~ qq state with the f /2~ gluino gives a 1/2+ i?-meson. The 
combination of the f + qq state with the gluino gives both a 1/2" nonet and a 3/2 + nonet. 
It is sufficient to focus only on the lightest states so we simplify the meson sector to include 
the two lightest nonets with J = 1/2. In [19] it is shown that as rrig — > oo, the mass splitting 
between states becomes constant: m^ — nip = 130 MeV, to# — nip = 40 MeV. It is important 
that the i?-pions (l/2 + ) are not significantly lower in mass than other non-strange i?-mesons, 
unlike the ordinary pions, because they are not pseudo-Goldstone bosons. Thus we take the 
i?-meson masses to be split only by the strange quark conteni[f|. With two light nonets and 
a gluino-ball we have 19 states: 10 degenerate non-strange i?-mesons, a gluino-ball we take 
to be degenerate, and a heavier set of 8 i?-kaons. We take the mass splitting between these 
states to be 130 MeV. For ease of reference we refer to all non-strange i?-mesons as tt from 
here on, and all strange i?-mesons as K. 

mi> - m~~ 130 MeV 



Table 1: The mass splitting of the i?-meson states stable against strong decays from [19]. 

In the i?-baryon sector the gluino combines with a qqq colour octet. The qqq state can have 
either J = l/2orJ = 3/2; combined with the spin- 1/2 gluino this gives a set of states with 
J — 0, 1, 2. Fermi statistics requires that the J = 3/2 qqq state be a flavour octet, while the 
J = 1/2 state can be a flavour singlet, octet or decuplet. Buccella et al [20] showed that for 
gluino masses above 1 GeV, the only i?-baryon states that are stable against strong decays 
are the J = flavour singlet guds (S°), the flavour octet with J = and the lighter of the 
two flavour octets with J = 1. 

We can estimate the weak-decay lifetime for Rn — >• S° + 7r ± to be 

^«1.6xlO-sec( I ^| v ) 2 , (1) 

by scaling the rate for A — > N + tt with particle masses as appropriate for 2-body decay 
kinematics. Thus, we can ignore weak decays in the rest of our analysis. 

The splittings between the states are summarised in Table [ffil. In practice we can ignore many 

3 The ft-pions are actually expected to be 40 MeV heavier than the ij-rho according to |19| . but this is 
an unnecessary level of refinement and we ignore it here. 

4 There is some uncertainty in the mass splitting that arises from the value of the zero point energy taken 
in the bag model calculation. We take the values in Table [2] for our model and then vary the splitting to 
check that our results are robust with respect to the splittings. 
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- m lo 


250 MeV 


m 8 


- m 8l 


130 MeV 



Table 2: The mass splitting of the i?-baryon states stable against strong decays from |20) . 

of the i?-baryons. Any singly strange state in the octet can decay strongly to the flavour 
singlet and so should not be considered to be quasi-stable. Moreover because the i?-baryons 
are created through interactions of the form tt + N — > N + ir we can ignore the doubly strange 
states: there is no way of creating a doubly strange state from the interactions of an i?-meson 
with the nucleons in the detector. We designate the singlet state as S°, the spin-0 states as 
h, p and the spin-1 states as A , A + . 

The mass calculations of [TH] and [2UJ were performed using several choices of bag model 
parameters. The primary effect of the different parameter choices is to alter the overall 
scale of the masses. Therefore we keep the interstate splitting within the meson and baryon 
sectors. We take the mass splitting between the meson and baryon sectors to be equivalent 
to the mass of the extra valence quark which we set to 330 MeV. In our final results we allow 
for this to be varied and find that this splitting has little effect on the overall conclusions. 



i?-baryons 


Ruds '■ 


m- a 


+ 280 MeV 


Rs-o: A+/° 


m- g 


+ 530 MeV 


i? 8 _ : h, p 


m g 


+ 660 MeV 


i?-mesons 




m-g 


+ 330 MeV 


k 


m g 


+ 460 MeV 



Table 3: Masses used for the i?-hadrons. 

Finally, we must set the overall mass normalisation with respect to the gluino mass. This 
value is sensitive to the details of the bag model calculations. We follow Buccella et al [20] , 
setting m(n) — m(g) + 660 MeV. This fixes the overall mass scale and gives the spectrum 
listed in Table [3l In our final results we allow this overall scale to vary and include the effect 
in our estimation of the theoretical uncertainty 

3 Propagation of R-hadrons through a detector 

Following production in a hard collision, a gluino would fragment dominantly into either an 
i?-meson or gluino-gluon state. This particle would then propagate through the detector. In 
order of increasing distance from the interaction point, a typical multipurpose experiment 
at a collider comprises an inner tracking system, a calorimeter and an outer set of muon 
chambers. If charged, the momentum of the i?-hadron would be measured in the inner and 
outer tracking systems. Hadronic and electromagnetic energy loss of the i?-hadron would be 
recorded in the calorimeter. 

Hadronic interactions of i?-hadrons occur predominantly due to the interactions between 
target nucleons and the light quark system accompanying the gluino, with the much more 
massive gluino acting as a spectator [TB]. The full list of hadronic processes we consider 
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is given in Appendix [X] For i?-hadrons produced at the Tevatron the light quark system 
typically has O(GeV) kinetic energy, implying that the interactions with nucleons should 
resemble low energy hadronic interactions. In such processes quarks can be exchanged with 
the quarks in the target nucleon, which can change the charge of the i?-hadron. Another 
process in i?-hadron scattering with nucleons is the conversion of i?-meson and gluino-gluon 
states to baryon states with pion production. Owing to lightness of the pions, such processes 
are exothermic and thus energetically favoured. Furthermore, the absence of pions in the 
scattering material renders the reverse process unlikely. In the scattering model used in this 
work, an _R-hadron undergoes on average ~ 5 hadronic interactions in 2 m of iron |13) , which 
is representative of the material within a typical calorimeter at a collider experiment. Thus 
most i?-hadrons would enter the muon detector as baryons, predominantly the spin-1 octet 
states R U ud and R u dd and the spin- zero singlet state S° (R u ds), as discussed in Sec. [2] 



4 Simulation of R-hadrons 



To simulate the iZ-hadron signal in the detector we employ several computer codes. 
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Figure 1: NLO gluino production cross sections at the tevatron for five different gluino 
masses. The dotted line shows the cross-section for direct production of a pair of gluinos. 
The solid line shows the total cross-section for the production of any pair of coloured SUSY 
particles, relevant because squarks decay virtually instantaneously to a gluino. As the squark 
masses increase the squarks decouple and the two cross-sections converge. 

Gluino production 

The cross sections for gluino pair production are calculated using Prospino2.1 |25j assuming 
all squarks have a common mass. The cross section dependence on the squark and gluino 
masses is shown in Figure [1] for the Tevatron. Two cross-sections are plotted for each gluino 
mass. The dotted line shows the direct pair production cross section for two gluinos. The 
solid line shows the total cross-section for the production of any pair of coloured SUSY states 
(squark-antisquark, gluino-(anti)squark, and gluino- gluino). As the squark mass increases, 
squark production falls off and the two cross-sections converge. Any squark produced decays 
to a quark and gluino, effectively instantaneously. We use this to defines three regimes in 

5 We ignore such small mass splittings, Am <C 1 GeV, that the squark decay is not effectively instantaneous; 
even though with small mass splitting the stop, sbottom and scharm states would be stable, these provide a 
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which we set mass limits: 

1. Decoupled squarks. For large squark masses the cross-section for the production of two 
gluinos becomes independent of the squark mass and is just the cross-section for the 
direct production of a gluino pair. This can be compared directly with the split SUSY 
scenarios commonly used in gluino searches. 

2. Light squarks. In the other extreme, that squarks are light (but heavier than the gluino 
as required for a long-lived gluino), an important source of gluino production is the 
production and decay of squarks, particularly through the process qq — > qq mediated 
by t-channel gluino exchange. The (anti-) squarks decay essentially instantly to a gluino 
and an (anti-) quarljf]. In this limit the (anti-) quark is soft and does not create a jet 
that could mask the R-hadron track and the gluino has the same kinematics as its 
parent squark. In this regime, all coloured SUSY particles can simply be replaced by 
a gluino for the purpose of this analysis and we take the cross-section for gluino pair 
production to be the summed cross-section to any pair of SUSY coloured states. 

3. Squark mass agnostic. When squarks are sufficiently heavier than the gluino, the quark 
emitted in the squark decay forms a jet which may obscure detection of the gluino 
track, thus reducing the sensitivity to gluinos produced via squark decay. The most 
conservative limit for an arbitrary squark mass, which can be given without detailed 
consideration of event properties, is therefore obtained by considering only directly 
produced gluino pairs. We simplify the parameter space by taking the light squark 
masses to be degenerat^, then use the squark mass value that minimizes gluino pair 
production; this occurs when the squarks are approximately degenerate with the gluino. 

Gluino hadronization 

Monte Carlo samples of the final states were generated for a y/s = 1.96 TeV pp collider 
using Pythia [15] ■ The squarks and gluinos are hadronized using the same hadronization 
scheme as in previous studies [321 [17] . Hadronization is restricted to R- mesons only, since R- 
baryon production should be suppressed in a similar way that ordinary baryon production is 
suppressed compared to meson production. As the proportion of baryons produced compared 
to all particles is small in data (< 5%), and we expect a similar ratio for i?-baryons, this will 
be counted as a systematic error on the cross section. The gluinoball formation probability 
is set to 10% compared to the i?-mesons. 

i?-hadron propagation 

This sample of gluinoballs and i?-mesons is then injected and propagated through the detec- 
tor. For this work, a model of i?-hadron scattering implemented in Geant-4 [22] is used [P7] . 
In view of the inherent uncertainties associated with modeling i?-hadron scattering, a prag- 
matic approach based on analogy with observed low energy hadron scattering is adopted. 
The scattering rate is estimated using a constant geometric cross-section of 12 mb per light 
(u, d) quark and 6 mb per strange quark, all 2-to-2 and 2-to-3 processes are allowed if they 
are kinematically feasible, and charge conservation is respected. The proportion of 2-to-2 
and 2-to-3 reactions is governed by phase space factors; no explicit constraints are applied 
to the probability of baryon number exchange. 

negligible contribution to the final states. 

6 We assume Am = 1 GeV for definiteness. There is no symmetry enforcing rrig = m sq , so consideration 
of smaller mass splittings (which require careful treatment of non-perturbative QCD) is not motivated. 

7 To an excellent approximation, the gluino pair production cross-section depends on the mass of the up 
and down squarks only. The up and down squarks are typically degenerate at the electroweak scale, so the 
assumption of a common squark mass is well founded. 
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Figure 2: Flavour fractions for i?-hadrons as a function of penetration depth in iron. The 
fractions are calculated from the percentage of gluino i?-hadrons that have not stopped at 
the given depth. The stopping fraction is shown by the long dashed (red) line. We label the 
populations that are significant past 2m to aid in differentiating these lines. 



Figure [2] shows the evolution of the flavour composition of an i?-hadron sample for a gluino 
mass of 200 GeV. A typical gluino hadronises to a meson immediately after production, 
converts to an R-baryon within the first meter or so of the detector, and then propagates 
through the calorimeter losing energy through electromagnetic interactions and changing 
charge through strong interactions with the valence quarks. Some fraction of the R-baryons 
stop in the detector, some enter the outer muon chamber charged, and the majority enter the 
muon chamber neutral. The kinematic distribution of the i?-hadrons is taken from a Pythia 
event sample for this gluino mass. 



5 Re-evaluating the CDF CHAMP limit for the case of 
meta-stable gluinos 

In this section we apply the i?-hadron production and propagation models developed above 
to obtain a limit on the mass of a quasi-stable gluino from the CDF Charge Massive Particle 
(CHAMP) search [BJ. The CDF detector is described in detail in [13]. For the purpose of 
this paper, we will be focusing on the number of nuclear interaction lengths between the 
interaction point and the muon system. This is depicted in Figure |31 The pseudorapidity 
region considered in [BJ is bounded by the grey dashed line. 

To estimate the response of the CDF detector, the model described in Sec. |3Jis supplemented 
with a sample of gluino i?-hadrons generated with Pythia |15j . Kinematic distributions were 
derived from the sample, and high-statistics simulation samples were generated obeying kine- 
matic distributions extracted from Pythia. The CDF detector was simulated using a depth 
of iron extracted from Figure [3] at the generated pseudorapidity. The flavour distribution of 
the emerging i?-hadrons is shown in Figure |3J 

To emulate the analysis performed in [BJ, only i?-hadrons 
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Figure 3: Matter distribution in the CDF detector as a function of pseudo-rapidity \rj\. 



10" 2 



1Q- 3 



- Stoppage: 5.4% 




: 


ii 



K* ft" K - SO P 1 



Figure 4: Flavour fractions for i?-hadrons emerging from the simulated geometry. 

• in the pseudo-rapidity interval |?y| < 0.7 

• in the (3 interval 0.4 < j3 < 0.9 

• with a transverse momentum above 40 GeV after traversal of the simulated iron 

were considered. Furthermore only i?-hadrons that were positively charged both immedi- 
ately after hadronization and after traversing the simulated geometry were considered. This 
satisfies the requirement in [6] that the R-hadron leaves a charged track in the central tracker 
and the outer muon chambers. 

The probability for an i?-hadron to be accepted turns out to be of order 10%. Denoting the 
single-object acceptance probability by x, the resulting efficiency for pair-production events 
to be accepted becomes 

e V is = 2x(l -x) +x 2 (2) 

yielding typical efficiencies of roughly 20%. eyis is dependent upon details of the model for 
production, hadronisation and scattering. To check that our results are robust with respect 
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Table 4: eyis for gluino pair-production events calculated at four mass points. Rows 1-6 
represent variations to the base model and are described in the text. 

to the theoretical uncertainty in these processes, we varied details of the model and studied 
the resulting change in eyis- Table |4] shows eyis calculated within the simulations for gluinos 
at four mass points. The model variations listed in the table are: 

1. Increasing all i?-hadron masses by 1 GeV. 

2. Increasing the i?-meson/i?-baryon mass splitting by 500 MeV. 

3. Setting the i?-meson/i?-baryon mass splitting to GeV. 

4. Multiplying the nuclear scattering cross section by a factor of 2 

5. Multiplying the nuclear scattering cross section by a factor of 0.5 

6. Using the Regge model from ref. fT3] to model the nuclear scattering. 

In 1-3 we explore the sensitivity of the results to uncertainty in the spectrum of H-hadron 
masses, while 4-6 explore the sensitivity to uncertainty in the hadronic interactions of R- 
hadrons. The only variation that results in a lower visible fraction arises from doubling the 
hadronic scattering cross-section. All other model variations result in a larger visible fraction 
that would give a stronger limit on the mass of the gluino than the base model we take 
here. Thus the limits we obtain with this analysis are both conservative and robust against 
theoretical uncertainties 

The result of the simulation is shown in Figure [5J The nearly horizontal line gives CDF's 
95%CL upper limit on stop squark pair production multiplied by the acceptance from [6]. 
The dependence of the acceptance on stop mass is virtually negligible, so can conservatively 
be extrapolated to the higher masses relevant for our limit as a constant. The red line 
represents a constant extrapolation. 

The three lines correspond to the three regimes outlined in section 2J In each case the 
cross-section is multiplied by the efficiency factors: 

8 Our results are also robust with respect to the ordering of the R-baryon mass hierarchy. The R-baryons 
entering the muon chambers are 23% charged to 77% neutral with the default spectrum (see Fig. due to 
the large fraction of R-baryons that convert to the light stable S°. If the S° were heavier, the R-baryons would 
be a democratic mixture of the lightest degenerate states - giving approximately 50% charged to 50% neutral 
when they enter the muon chambers. A factor 2.2 increase in the charged fraction doubles the predicted 
signal improving the final limits to 423 GeV, 430 GeV, 345 GeV for the decoupled, degenerate and agnostic 
squark mass scenarios respectively. Therefore the light S° is both well motivated, and provides a conservative 
limit. 
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Figure 5: Calculated gluino pair production cross sections with efficiencies applied compared 
to the Tevatron pair production limit (horizontal dashed line). The solid line, ctes, corre- 
sponds to the split SUSY limit with decoupled squarks. The dotted line, cr' ES , corresponds 
to the squark mass agnostic scenario. The dot-dashed line, cr^g, corresponds to the scenario 
with squarks degenerate with the gluino, in which the gluinos are primarily formed through 
the decays of (anti-)squarks. 

• Acceptance eyis-' The application of the nominal acceptance factors from Table [4] 

• Reconstruction efficiency £r C c' A blanket reconstruction efficiency of 38 % is applied 
in accordance with [6] . 

The solid line, (Jes, corresponds to the split SUSY limit with decoupled squarks. The dotted 
line, <r Eff , corresponds to the squark mass agnostic scenario. The dot-dashed line, <j'^ s , 
corresponds to the scenario with squarks degenerate with the gluino, in which the gluinos 
are primarily formed through the decays of (anti-) squarks. 

The grey bands represent the systematic uncertainties. To calculate the error the renormal- 
ization and factorization scales were increased and decreased by a factor of 2 and the standard 
CTEQ6 pdf [25] was replaced by the MSTW pdf [27]. The uncertainty is dominated largely 
by the scale variation, and the error band represents the sum in quadrature of the two sources 
of error. This pQCD systematic uncertainty far exceeds the uncertainties quoted in Table |4] 
and the latter are omitted in Figured] 

We take the mass where the 95 %CL Tevatron limit crosses the lower error band of the 
effective cross section <7Eff to be a lower limit on the mass of a stable gluino. This gives lower 
limit on the gluino mass of 397 GeV for decoupled squarks, 407 GeV for squarks degenerate 
with the gluino, and 322 GeV in the squark mass agnostic case. 

6 Conclusions 

We have used the results of the search for stop squark i?-hadrons by CDF [6] to derive a 
limit on the production of a gluino-containing i?-hadron which traverses the detector without 
decaying, and hence a lower limit on the mass of a gluino with r > 30 ns. Starting from 
QCD predictions for the i?-hadron mass spectra, we model the production, hadronisation and 
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scattering of i?-hadrons within the CDF detector and calculate the fraction of i?-hadrons that 
mimic the signal investigated in [6] . We estimate the effects of the uncertainties in the mass 
spectrum, gluino production cross section and details of the scattering model. After including 
systematic uncertainties, we find a lower limit of 407 GeV on the mass of a quasi-stable gluino 
whose mass is approximately degenerate with squarks, a limit of 397 GeV on the gluino mass 
if squarks are much more massive than gluinos, and a limit of 322 GeV that is valid for 
any (degenerate) squark mass; gluino mass limits can readily be derived for other squark 
mass scenarios. These are currently the strongest experimental limits on the mass of a stable 
or long-lived gluino and show the continued strength of Tevatron data in constraining new 
physics. This work also lays a foundation for interpreting future LHC searches. 
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A List of hadronic processes 

This appendix contains the full process lists used in this paper. No distinction is used 
between i?§-o and i?s-i -R-baryon states in this list, excluding processes that specifically 
mix the states. It is understood that any process written down for i?s~o particles is also 
valid for i?s-i- 
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# 



uds 



N -> R uds + N 



S° 



• S u + n 
S°+p 



# 8 - + N -> flg-i + JV 



# 



uds 



N -> R 



uds 



N 



5^ 

5° 

5° 



n 
n 
P - 



TT 



- 5"' + n + tt 
S°+ p + tt° 



S° S° + n + TT- 



Ruds + N 

S v + n^ 
S°+n 
S° + n 
S°+p- 
S°+ P - 
S°+ P - 



Rs o + K + TV 

• n + if u + n 

• n + if - + p 

• p + K~ + n 
p + K° + n 
p + K~ +p 
h + K° +p 



#8 



+ N 



R 



18-0 



+ N 



h+p 
h+p 

p + n- 
p + n - 
p + p- 



> n + n 
•n+p 
•p + n 
•p + n 
■h+p 
p + p 



R 



8-0 



n + n 
h + n 
h + n 
h+p- 
h+p- 
h+p- 
h+p- 
p + n- 
p + n - 
p + n - 
p + n- 
p + p- 
p + p- 
p + p- 



N -» fis-o + N + tt 

■> n + n + 7T U 

■> n + p + 7r~ 

■> p + n + tt~ 
o 



h + p + TT 
n + n + 7r 4 
p + n + ir' 
p + p + ir~ 
p + n + tt' 



o 



o 



n 

h - 



P- 
■ n ■ 



- p + TT 

- n + tt 

-p + TT 







• p" + p + 7T~ 

n 



h+p 
h+p 
p + n 
p + n 
P + P 



#8-0 + N 



A + n 
A°+p 
A+ + n 
A+ + n 
A°+p 
A+ +p 



R 8 -!+N + n 

.0 



h + n^~AP +n + it u 
h + n^A°+p + ir~ 
n — » A+ + n + 7T 
A°+p + 7r° 
A + n + 7T+ 
A+ + n + tt 
A+ + p + 7^- 
A + + n + tt° 

.0 



n 

h +p -> A 

n+p 

P 
n 
n 
n 



p + n 
p + p- 
p + p 



#8-0 + W 



► A° + p + 7T° 

► A° + n + 7T+ 

► A+ + p + 7T" 

- A+ + p + 7T° 

- A+ + n + 7T+ 

A° + P + 7T + 



n + n 
n+p 
n+p 
p + n 
p + n 
p + p- 



# M ds + K + N 
K u + n 



> S° + K + p 

> S 10 + if + + n 

> 5° + #°+p 

■ 5° + K+ + n 
- S° + K++p 



Table 5: Process lists galore for baryons 
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N -> 7T + N 



- n — >• 7r u + n 



+ n — »■ 7r~ + p 

+ p -» 7T° + p 

+ p — » 7T+ + n 
" + n — > 7r~ + n 

" + P -> 7T~ + P 

" + p -> 7r° + n 



+ n — > 7r + + n 
+ n — > 7r° + p 

+ p — > 7T+ + p 



A/" 



7T + TV + 7T 



"IF 

7^° 



- n —> 7i- u + n + 7T U 



n 
n 



7T° + P + 7T 



■ Tl — > 7T + p + 7T l; 



7f° 



7f° 



+ n 
+ P- 
+ P- 
+ P- 
+ P- 
+ P- 



7T° + P + 7T 
7T° + Tl + 7T 4 
7T+ + n + 7T 1 
7T+ + p + 7T~ 
7T~ + p + 7T+ 







7T° 



7T 

7T" 

7T~ 



- AT fig-Q + 7T 

+ n — >• n + 7r u 
+ n — > p + 7r~ 
+ p — > n + 7r + 

+ p -> p + 7T° 

" + n -> p + 7T° 
" + n — > n + 7T+ 

" +p — > p + 7T + 

" + n — > n + 7r~ 
" + p — > n + 7r° 

" +p -> p + 7T~ 



7T + AT -)• AVo +TT + TT 



7T U + n 
7T° + n 
7T° + n 
fr + p- 
fr° +p - 

7T° + p - 

7r + + n 
7T+ + n 
7r + + n 

7T + + p 
7T + + p 

7r~ + n 
7r~ + n 

7T~ +P 
7T~ + p 
7T~ +p 



n 
n 



p + 7T + 7T 



->• U + TT + 71" 
-)• p + 7T° + 7T° 
-> p + 7T + + 7T~ 
-)• p + 7T° + 7T° 
— > p + 7T+ + 7T~ 
-)• U + 7T+ + 7T° 
— > p + 7T + + 7T~ 
— > n + 7T+ + 7T+ 
— > n + 7T~ + 7T 
— » p + 7T~ + 7T~ 

-> fl + 7T° + 7T° 

+ 







n + 7r 1 



■ 7T 



— > p + 7T + 7T 







7T + N 



7r u + n 

7T° + p - 

7T+ + n 

7T~ + p 



7T + N 



7f u + n - 
7T° + n - 
n° +p - 

7f° +p - 

7r + + n - 
7T+ + n 

7T+ + p - 

7T _ + n ■ 

7T~ + p - 
7T _ + P - 



Ruds 
^1 



+ K- 
^4 



- 7T 

S 10 + + 7T- 

5° + tf + + 7T° 
5° + K° + TT+ 
S"+K++ TT° 
S"+K°+ TT+ 
S°+K+ + 7T+ 
S° + K° + 7T- 

S° + K + + 7T- 



K + N ^ K + N 



K+ + 
K+ + n ^ 
K+ + p ^ 
r+m 
+p -> 
X" +p -> 
K a + 
K a +p^> 
K°+p^ 
K° + n^ 
K° + n^ 
K°+p^ 



K+ +n 

k° + P 

K++p 

k- 

K~ 

k° 
k° 

k a - 
k+ 

K a + n 
K-+p 

k° + P 



+ n 

- n 
n 
P 

- n 



K + N^K + N + tt 



K + +n^K + +n + i\ 
k + + n K + p + V 
X++n^X n + n + 7r 



n-0 



^0 



/?+ + n -)• + p + TT" 

+ p -)• + p + 7T° 

/?+ + p ->• + n + 7T+ 

+ p -)• K° + p + 7T+ 

K~ + n -)• A" - + n + 7T° 
X" + n ->• A" - + p + 7r~ 
X" + n ->• AT + n + 
A > - + p -)• AT" + p + 7T° 
A > " + p -)• K~ + n + 7T+ 
k-+p^k° + n + TT 

k-+p^k°+p + TT 

K Q + n -)• X + n + TT' 
AT + n -4 X + p + 7T- 
K a + n ^ K+ + n + 

k° + P -+ k° + P + it" 

K Q K° +H + TT+ 

k°+p^k + +n + 7T 
k°+p^k + +p + n 
K + n^k° + n + TT 
K° + n -> K° + p + n- 
K° + n -> K- + p + tt° 
K° + n ->• X" + n -J- 
X + p X + p + n 
A; + p -4 X + n 

k° + p^- k~ +p + tt + 



™-0 



- 7T 
T° 



Table 6: Process lists galore for mesons 
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K + N -> 


Ruds + 7T 




if - +n - 


>S° + 7T- 




K~ +p- 


>S° + 7T° 




K° + n^ 


5° + 7T 




k° + P ^ 


S°+7T+ 




K + N ^ R uds + 7T H 


- 7T 


if" +n - 


► 5° + 7T- - 




A"" + p -> S*° + 7T° + 7T° 


A"" + p -> 5° + 7T+ H 


- 7T 




5° + 7T° + 


7T° 


K° + n^ 


5° + 7T+ 4 




k° + P ^ 


S° + 7T+ + 


7T 


K + N -> 


#8-0 + A" 




K+ +n- 


> n + if+ 




K+ +n- 


>p + K° 




k+ + p — 


■p + K+ 




k+ + p — 


>n + if° 




K~ + n- 


>h + K~ 




R- +p- 


>p + K~ 




k- +p- 


>n + K° 




K° + n^ 


n + K° 




k° + p^ 


p + K° 




k° + p^ 


n + K+ 




K° + n^ 


n + K° 




K° + n^ 


p + K- 




k° + p^ 


P + k° 





Table 7: Baryon number changi 



K + N -» i^p + X + 7T 

if + + n -> n + if + + 7T° 
if + + n -)■ n + K° + 7T+ 
if + + n -> p + K° + tt° 
K+ + n -> p + if + + 7T" 

if + + p -)• p + if + + 7T° 
if + + p -> p + if° + 7T+ 
if + + p -)• U + if + + 7T+ 



if- 


+ n 


-+h + K~ 


+ 7T 


if- 


+ n 


-^h + k° 


+ 7T- 


R- 


+ n 


-^p + K- 


+ 7T- 


k~ 


+ P 


-^p + K~ 


+ 7T° 


k~ 


+ P 


-^p + k° 


+ 7T- 


R- 


+ P 


-> n + if + tt° 


R- 


+ P 


-> n + if" 


+ 7T+ 


k° 


+ n 


^ n + if° 


+ 7T° 


k° 


+ n 


-> n + if- 


+ 7T+ 


k° 


+ n -> p + if 


+ 7T° 


k° 


+ n 


^p + if° 


+ 7T" 


k° 


+ P- 


-^p + k 


+ 7T° 


k° 


+ P- 


p + K~ 


+ 7T+ 


k° 


+ P- 


-> n + if° 


+ 7T+ 


k° 


+ n 


-)• ra + if" 


+ 7T° 


k° 


+ n 


^ n + if° 


+ 7T" 


k° 


+ n 


^p + if- 


+ 7T- 


k° 


+ P- 


^p + if- 


+ 7T° 


k° 


+ P- 


-^p + if° 


-t-TT" 


k° 


+ P- 


-> n + if° 


+ 7T 


k° 


+ P- 


-> n + if" 


+ 7T+ 



processes for strange mesons 
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LHC Vs=7 TeV 
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